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ABSTRACT 

We report the discovery by the CoRoT satellite of a new transiting giant planet in a 2.83 days orbit about a V- 
1.08 ± 0.08 M ,R, = 1.1 ±0.1 R , T eB = 5675 ± 80 K). This new planet, CoRoT- 12b, has amass of 0.92: 
Rj U p. Its low density can be explained by standard models for irradiated planets. 
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1. Introduction 

Because of their special geometric configuration, a wealth of in- 
formation can be learned about transiting extrasolar planets (e.g., 
Winn 2010), making them very important for our understand- 
ing of the vast planetary population hosted by our galaxy. They 
are the only exoplanets for which accurate measurements of the 
mass and radius are available. Furthermore, their atmospheric 
properties can be studied during their transits and occultations 
(e.g., Deming & Seager 2009). 

More than 70 extrasolar planets transiting their parent 
stars are now knownj, most of which having been discovered 
by dedicated photometric surveys. Among these, the CoRoT 
(Convection, Rotation, and planetary Transits) space mission 
(Baglin et al. 2009) stands out as a pionner project. Because 
of its excellent instrumental capabilities and its low Earth or- 
bit, CoRoT can monitor the same fields of view with a very high 
photometric precision for up to five months. This makes possible 
the detection of planets that would be out of reach for ground- 
based surveys, as demonstrated for instance by its discovery of 
the first transiting 'Super-Earth' CoRoT-7b (Leger et al. 2009; 
Queloz et al. 2009), and the first 'temperate' transiting gaseous 
planet CoRoT-9b (Deeg et al. 2010). 

We report here the discovery of a new planet by CoRoT, a 
'hot Jupiter' called CoRoT-12b that transits a my - 15.5 so- 
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* The CoRoT space mission, launched on December 27, 2006, has 
been developed and is operated by CNES, with the contribution of 
Austria, Belgium, Brazil , ES A (RSSD and Science Program), Germany 
and Spain. 

1 See, e.g., Jean Schneider's Extrasolar Planet Encyclopedia at 
http://exoplanet.eu 



lar analog star. We present the CoRoT discovery photometry in 
Sec. 2. The follow-up, ground-based observations establishing 
the planetary nature of CoRoT-12b are presented in Sec. 3, while 
the spectroscopic determination of the parameters of the host star 
is described in Sec. 4. A global Bayesian analysis of the CoRoT 
and follow-up data is presented with its results in Sec. 5. Finally, 
we discuss the inferred properties of the CoRoT-12 system in 
Sec. 6. 



2. CoRoT photometric observations 

Table 1 presents the ID, coordinates and magnitude of CoRoT- 
12. This star is located in a field near the galactic anti-center 
direction, in the Monoceros constellation. It was monitored by 
CoRoT from October 24, 2007 to March 3, 2008 (CoRoT run 
LRaOl; see Rauer et al. 2009, Carone et al. in prep.). 

The transits of CoRoT-12b were noticed after 29 days by the 
so-called 'alarm mode' pipeline (Surace et al. 2008). The time- 
sampling was then changed from 512s, the nominal value, to 
32s. The processed light curve (LC) of CoRoT-12 is shown in 
Fig. 1 . This monochromatic LC consists of 258 043 photometric 
measurements for a total duration of 131 days. It results from 
the processing of the raw CoRoT measurements by the stan- 
dard CoRoT pipeline (version 2.1, see Auvergne et al. 2009), 
followed by a further processing (outliers rejection and system- 
atics correction) similar to what is described by, e.g., Barge et 
al. (2008) and Alonso et al. (2008). 47 transits of CoRoT-12b 
are present in the LC, 36 of them being found in its oversam- 
pled part. Some discontinuities are present in the LC. They were 
caused by energetic particles hits during the crossings of the 
South-Atlantic Anomaly by the satellite. A large jump of the 
measured flux (more than 5%) caused by the impact of a cosmic 
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CoRoT window ID 
CoRoT ID 
UCAC2 ID 
USNO-A2 ID 
USNO-B1ID 
2MASS ID 
GSC2.3 ID 

Coordinates 



LRaOl E2 3459 

0102671819 

31290403 

0825-03015398 

0887-0101512 

J06430476-01 17471 

SB3BK006251 



RA (J2000) 
Dec (J2000) 

Magnitudes 



06 43 03.76 
-01 17 47.12 



Filter 



Mag 



Error 



B" 

V" 
r' a 
i' a 

J" 
H" 

K" 



16.343 
15.515 
15.211 
14.685 
14.024 
13.630 
13.557 



0.080 
0.052 
0.040 
0.069 
0.029 
0.030 
0.041 



" Provided by Exo-Dat (Deleuil et al, 2009). 
b from 2MASS catalog (Skrustkie et al. 2006). 

Table 1. IDs, coordinates and magnitudes for the star CoRoT- 
12. 



ray on the detector can also be noticed in the last part of the LC. 
The processed LC shown in Fig. 1 has an excellent duty cycle of 
91%. 

Despite that its CoRoT LC shows some kind of irregular 
variations with a peak-to-peak amplitude of 2.3%, CoRoT-12 
appears to be a photometrically quiet star. Except for the tran- 
sit signal (see below), the discrete Fourier-transform of the LC 
shows no clear periodicity over the noise level. The rotational 
period of the star cannot thus be constrained from the CoRoT 
photometry. 

Periodic transit-like signals are clearly visible in the LC, as 
can be seen in Fig. 1. Initial values for the orbital period P 
and transit epoch Tq were determined by trapezoidal fitting of 
the transit centers, as described by Alonso et al. (2008). The 
resulting values were T = 2545398.6305 ± 0.0002 HJD and 
P - 2.82805 + 0.00005 days. These values were used to sched- 
ule the ground-based follow-up observations (see next Section), 
and also as initial guesses for the global analysis presented in 
Sec. 5. 



3. Ground-based observations 

The following ground-based observations were performed to es- 
tablish the planetary nature of CoRoT-12b and to better charac- 
terize the system. 
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Fig. 1. Top: Normalized CoRoT LC of the star CoRoT-12. The 
oversampled part of the LC was binned to the same time bin than 
its first part for the sake of clarity. Bottom: zoom on a transit of 
CoRoT-12b. 
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Fig. 2. The sky area around CoRoT-12 (brightest star near the 
centre). Left: R-fi\Xer image with a resolution of 1.3" taken with 
the INT/WFC. Right: Image taken by CoRoT, at the same scale 
and orientation. The jagged outline in its center is the photo- 
metric aperture mask; indicated are also CoRoT's x and y im- 
age coordinates and positions of nearby stars from the Exo-Dat 
(Deleuil et al 2009) database. 



3.1. Imaging - contamination 

CoRoT has a rather poor optical resolution, so performing high- 
resolution ground-based imaging of its fields is important, not 
only to assess the possibility that the eclipse signals detected by 
CoRoT are due to contaminating eclipsing binaries, but also to 
estimate the dilution of the eclipses measured by CoRoT caused 
by contaminating stars (see Deeg et al. 2009 for details). 

Imaging of the target field was undertaken with the 2.5m INT 
telescope during pre-launch preparations (Deleuil et al., 2009) 
and with the IAC80 telescope during candidate follow-up (Deeg 



et al. 2009). It found no nearby contaminating star that could be 
a potential false alarm source, i.e. that mimiks CoRoT's signal 
while being an eclipsing binary star (see Fig. 2). 

Using the method describe by Deeg et al. (2009), the fraction 
of contamination in the CoRoT-12 photometric aperture mask 
was estimated to be 3.3 + 0.5%. It is mostly due to a 3.5 mag 
fainter star that is 8.5" SW. This small dilution was taken into 
account in our analysis presented in Sec. 5. 
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3.2. Radial velocities - spectroscopy 

Four radial velocity (RV) measurements were obtained with the 
HARPS spectrograph (Pepe et al. 2002b, Mayor et al. 2003) on 
the 3.6-m telescope at ESO La Silla Observatory (Chile), on 
October 2008 (HARPS program 082.C-0120). These first data 
were made using the high efficiency mode EGGS in order to 
establish the planetary nature of the companion, showing a de- 
tectable and low-amplitude radial velocity variation in phase 
with the CoRoT ephemeris, with the shortest exposure time. Ten 
additionnal measurements were recorded with HARPS, from 
November 27, 2009 to February 05, 2010 (HARPS program 
184.C-0639). These newer data points were acquired using the 
high accuracy mode HAM to increase the precision of the RV 
measurements compared to the about 30 m s~' of systematic er- 
rors of the high efficiency mode (Moutou et al. 2009), and with- 
out simultaneous thorium (obj_AB mode) in order to monitor 
the Moon background light on the second fiber B. Radial ve- 
locities were obtained from the HARPS spectra by computing 
weighted cross-correlation with a numerical G2 mask (Baranne 
et al. 1996; Pepe et al. 2002a). 

Sixteen spectra of CoRoT-12 were also acquired with 
the HIRES spectrograph on the Keck I telescope as part of 
NASA's key science project in support of the CoRoT mission. 
Differential RVs were computed from these spectra with the 
Austral code (Endl et al. 2000). First, ten spectra were gathered 
during a transit of CoRoT-12b in January 2009. Unfortunately, 
the used set-up of the slit decker did not allow a proper subtrac- 
tion of the sky background, leading to RV systematics with an 
amplitude of a few dozens m s~' , so we decided to reject these 
data. Six other HIRES RVs were obtained between December 
2009 and January 2010. For these six spectra, the set-up of the 
slit decker was changed, leading to a proper background subtrac- 
tion. 

Our HARPS and HIRES measurements are presented in 
Table 2. An orbital analysis was performed treating the three 
sets of RV measurements (HARPS HAM, HARPS EGGS, and 
HIRES) as independent data sets with different zero point ve- 
locities. The orbital solution was made keeping the period and 
ephemeris fixed to the CoRoT values, but allowing the zero point 
offsets to be fit in a least square way. Fig. 3 shows the resulting 
orbital solution which is in phase with the CoRoT photomet- 
ric signal. The resulting eccentricity (0.03 ± 0.13) was consis- 
tent with zero while the semi-amplitude was 124+15 ms 4 . 
Assuming a solar-mass host star, this semi-amplitude translates 
into a transiting object with a mass of about 0.9 Mj. In Sec. 5 
we present a revised orbit obtained using a global analysis. 

The residual RVs were analyzed after removing the orbit to 
look for the possible presence of additional companions. No 
significant variations were found, but given the sparseness of 
the measurements we cannot exclude the presence of additional 
companions with a good confidence. 

The HARPS cross-correlation functions were analyzed using 
the line-bisector technique (Queloz et al. 2000). Fig. 4 shows the 
correlation between the bisector and RV measurements. The cor- 
relation coeffient of all the RV-bisector measurements, r, has a 
value of 0.56 with a probability 0.026 that the data is uncorre- 
cted. Ostensibly this correlation looks to be significant, but we 
do not believe that to be the case as this correlation is largely 
driven by one outlier in the HARPS data and another EGGS 
measurement. When one examines only the HARPS data the co- 
efficient drops to r — 0.47 with a probabilty of 0. 15 that the data 
is uncorrected. Removing one outlier point lowers the correla- 



tion coefficient to r — 0.32 with the probability of no correlation 
being 0.37. 

We believe that the modestly high correlation coefficient may 
be an artifact of the bisector error being more than a factor of two 
larger than the RV measurement error and the paucity of mea- 
surements. To test this we generated fake bisector/RV data con- 
sisting only of random noise that was sampled the same way as 
the real data. The standard deviations of the fake measurements 
were consistent with the median error of the RV and bisector 
measurements. In approximately 40% of the cases the correla- 
tion coefficient of these random data had correlation coefficients 
at least as large as that of the real data. The RV-bisector correla- 
tion coefficient that we measure is consistent with random noise 
coupled with sparse sampling. This discards the possibility that 
the periodic signal detected in these RVs is caused by a blended 
eclipsing binary. Taking into account the fact that CoRoT-12 is 
a solar analog star (see Sect. 4), we interpret thus the eclipses 
detected in CoRoT photometry as transits of a new giant planet, 
CoRoT-12b. 



HJD 


RV 


O-RV 


Bisector 


(days) 


(kms- 1 ) 


(kms-') 


(kms" 1 ) 


HARPS EGGS 


2454745.86036 


12.1740 


0.0221 


-0.0256 


2454746.83735 


11.9341 


0.0369 


-0.1307 


2454747.86641 


12.0904 


0.0198 


-0.0320 


2454763.81411 


11.9856 


0.0121 


-0.0435 


HARPS HAM 


2455163.73528 


12.1193 


0.0458 


0.0429 


2455165.71941 


11.9857 


0.0263 


-0.0310 


2455167.72180 


12.0570 


0.0342 


0.0105 


2455219.63940 


12.0051 


0.0195 


-0.0111 


2455220.68849 


12.2435 


0.0167 


-0.0188 


2455226.66329 


12.2126 


0.0168 


0.0246 


2455227.68971 


12.0348 


0.0442 


0.0335 


2455229.64150 


12.2355 


0.0292 


0.0259 


2455231.68894 


12.1640 


0.0240 


0.0167 


2455233.60091 


11.9993 


0.0292 


-0.0945 


HIRES 


2455170.99823 


-0.0573 


0.0286 




2455223.00984 


0.0490 


0.0144 




2455223.02060 


0.0466 


0.0144 




2455223.98643 


0.1633 


0.0143 




2455224.93395 


-0.0659 


0.0175 




2455224.94528 


-0.0979 


0.0220 





Table 2. HARPS and HIRES radial velocity measurements for 
CoRoT-12. The HARPS RVs are absolute, while the HIRES RVs 
are differential (measured relative to a stellar template). The bi- 
sectors were not measured from the HIRES spectra. 



4. Stellar parameters 

Two master spectra were used to determine the atmospheric pa- 
rameters of the star. The first of them was made by co-addition 
of the seven HARPS HAM spectra which were not strongly con- 
tamined by the Moon background light. The resulting master 
spectrum had a signal-to-noise ratio (SNR) about 40 in the con- 
tinuum. The second master spectrum was obtained from the co- 
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Fig. 3. HARPS and HIRES RVs phase-folded on the CoRoT 
ephemeris and overimposed on the best fit orbital model. 
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Fig. 4. Bisector versus RV measured from the HARPS spectra. 
Errors of twice the RV errors were adopted for all the bisector 
measurements. 



addition of two Keck spectra and had a SNR about 100 in the 
continuum. 

The methodology used to analyze these two master spectra 
was mainly based on the semi-automatic package VWA (Bruntt 
et al. 2002, 2008, 2010), and is thoroughly described by Deleuil 
et al. (2008) and Bruntt et al. (2010). The derived atmospheric 
parameters and elemental abundances are presented in Table 3. 

The Li I line at 670.78 nm was not detected in both mas- 
ter spectra, nor any hint of chromospheric activity. From this, 
the low rotational velocity measured in the spectra, and the low 
photometric variability noticed in the CoRoT LC, CoRoT-12 ap- 
pears thus to be a quiet and slowly rotating solar analog star. 

Using T e ff and log g from the VWA spectroscopic analy- 
sis, we estimated the absolute magnitude My — 4.75 mag and 
colour excess E(7 - K) ^ 0.08 mag from the Allen's tables (Cox 
2000). We calculated the corresponding interstellar absorption 
A v * 0.46 (using Ay = (5.82 ± 0.1) x E(7 - K); Cox 2000), to 
estimate, with the V apparent magnitude, the distance of the star 
tobed= 1150 + 85pc. 



Teg 


5675 ± 80 K 


l°gg 


4.52 ± 0.08 


^mic 


0.6 ±0.2 kms" 1 


^mac 


1.5 ±0.3 kms- 1 


v sin i 


1.0 ± 1.0 kms- 1 


d 


1150±85pc 


[Fe/H] 


0.16 ±0.10 


[Na/H] 


0.17 ±0.06 


[Mg/H] 


0.13 ±0.07 


[Al/H] 


0.15 ±0.10 


[Si/H] 


0.12 ±0.08 


[Ca/H] 


0.09 ±0.10 


[Sc/H] 


0.22 ±0.15 


[Ti/H] 


0.05 ± 0.09 


[V/H] 


0.02 ± 0.08 


[Cr/H] 


0.17 ±0.09 


[Mn/H] 


0.20 ±0.13 


[Co/H] 


0.16 ±0.14 


[Ni/H] 


0.21 ± 0.08 



Table 3. Stellar parameters and elemental abundances derived 
for CoRoT-12 from our VWA spectroscopic analysis. 



5. Global analysis 

5.1. Description 

We performed a thorough global analysis of the CoRoT transit 
photometry and HARPS/HIRES RVs to get the strongest con- 
straints on the system parameters. First, we cut the parts of the 
CoRoT LC located within 0.15 days of the transit mid-times 
deduced from the preliminary transit ephemeris presented in 
Sec. 2, getting thus 47 individual transit LCs. Considering their 
large number of measurements, we decided to stack the mea- 
surements of the 36 over-sampled transit LCs per 4, to speed-up 
our analysis. This binning did not affect our final precision on 
the system parameters, as the resulting folded LC (see Fig. 5) is 
still well sampled. 

Our analysis was done with the adaptative Markov Chain 
Monte-Carlo (MCMC) algorithm presented by Gillon et al. 
(2009; 2010). MCMC is a Bayesian inference method based 
on stochastic simulations that samples the posterior probabil- 
ity distributions of adjusted parameters for a given model. Our 
MCMC implementation uses the Metropolis-Hasting algorithm 
(see, e.g., Carlin & Louis 2008) to perform this sampling. Our 
nominal model was based on a star and a transiting planet on a 
Keplerian orbit about their center of mass. More specifically, we 
used a classical Keplerian model for the RVs, while we modeled 
the eclipse photometry with the photometric eclipse model of 
Mandel & Agol (2002) multiplied by a baseline model consist- 
ing of a different fourth-order time polynomial for each of the 
47 CoRoT time-series. The coefficients of these baseline models 
were determined by least-square minimization at each steps of 
the Markov chains (see Gillon et al. 2010 for details). 

Our analysis was composed of a nominal MCMC run, fol- 
lowed by two other MCMC runs having different specificities 
that are described below and summarized in Table 4. Each of 
the MCMC runs was composed of five Markov chains of 10 5 
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steps, the first 20 % of each chain being considered as its burn- 
in phase and discarded. For each run, the convergence of the five 
Markov chains was checked using the statistical test presented 
by German and Rubin (1992). 

The correlated noise present in the LCs was taken into ac- 
count as described by Gillon et al. (2010), i.e., a scaling factor 
was determined for each LC from the standard deviation of the 
binned and unbinned residuals of a preliminary MCMC analysis, 
and it was applied to the error bars (see also Winn et al. 2008). 
For the RVs, a 'jitter' noise of 5 m s was added quadratically 
to the error bars, this value being an upper limit for a quiet solar- 
type star like CoRoT-12 (Wright2005). Practically, this low jitter 
noise has no impact on the posterior distributions of the system 
parameters, as CoRoT-12 is faint and the RV precision is photon 
noise/background contamination limited. For the four HARPS 
measurements obtained with the EGGS mode, a systematic er- 
ror of 30 ms 4 was also added quadratically to the error bars 
(see Sec. 3.2). 

In all three MCMC runs, the following parameters were 
jump parameter^: the planet/star area ratio (R p /R s ) 2 , the transit 
width (from first to last contact) W, the parameter V — a cos i/R, 
(which is the transit impact parameter in case of a circular or- 
bit), the orbital period P and time of minimum light Tq, the 
two Lagrangian parameters e cos a> and e sin a> where e is the 
orbital eccentricity and u> is the argument of periastron, and the 
parameter K 2 = K Vl - e 2 P 1 ^, where K is the RV orbital semi- 
amplitude (see Gillon et al. 2009, 2010). We assumed a uniform 
prior distribution for all these jump parameters. To take into ac- 
count the small dilution of the signal due to contaminating stars 
(see Sec. 3.1), the jump parameters (R p /R s ) 2 was divided at each 
step of the MCMC by a number drawn from the distribution 
Af(1.033,0.005 2 ) before being used in the computation of the 
eclipse model. 

We did not assume a perfectly circular orbit in any of our 
MCMC runs despite that a circular orbit is compatible with the 
results of our orbital analysis of the RVs (see Sec. 3.2). Indeed, 
most short-period planets could keep a tiny but non-zero eccen- 
tricity during a major part of their lifetime (Jackson et al. 2008), 
so fixing the eccentricity to zero is not justified by tidal theory 
and could lead to overoptimistic error bars on the system param- 
eters. 

We assumed a quadratic limb-darkening law, and we allowed 
the quadratic coefficients u\ and u 2 to float in our MCMC runs, 
using as jump parameters not these coefficients themselves but 
the combinations c\ — 2 x u\ + u 2 and c 2 - u\ - 2 x u 2 to 
minimize the correlation of the obtained uncertainties (Holman 
et al. 2006). To obtain a limb-darkening solution consistent 
with theory, we decided to use normal prior distributions for u\ 
and u 2 based on theoretical values. Sing (2010) presented re- 
cently a grid of limb-darkening coefficients specially computed 
for the C0R0T non-standard bandpass and for several limb- 
darkening laws. We deduced the values u\ - 0.47 ± 0.03 and 
u 2 = 0.22 ± 0.02 from Sing's grid for the spectroscopic param- 
eters of CoRoT-12 and their errors (Table 3). The correspond- 
ing normal distributions N(0A7, 0.03 2 ) and N(0.22, 0.02 2 ) were 
used as prior distributions for u\ and u 2 in our MCMC analysis. 

At each step of the Markov chains, the stellar density de- 
duced from the jump parameters, and values for T e ff and [Fe/H] 
drawn from the normal distributions deduced from our spectro- 
scopic analysis, were used as input for the stellar mass calibra- 
tion law deduced by Torres et al. (2010) from well-constrained 



detached binary systemo Using the resulting stellar mass, the 
physical parameters of the system were then deduced from the 
jump parameters at each MCMC step. To account for the uncer- 
tainty on the parameters of the stellar calibration law, the values 
for these parameters were randomly drawn at each step of the 
Markov chains from the normal distribution presented by Torres 
et al. (2010). 

In our second MCMC run (labeled MCMC 2 in Table 4), we 
also used as data the parts of the C0R0T LC located within 0.2 
days of the occultation mid-times deduced from the best fit tran- 
sit ephemeris of our nominal MCMC run. The goal of this run 
was to obtain an upper limit for the depth of the occultation in 
the C0R0T photometry. For this run, the occultation depth was 
thus also a jump parameter. 

Finally, we assessed the perfect periodicity of the transits of 
CoRoT-12b in our third run (labeled MCMC 3 in Table 4). For 
this run, a transit timing variation (TTV) was considered as jump 
parameter for each of the 47 transits. Obviously, the orbital pe- 
riod could not be determined unambiguously without any prior 
on these TTVs, so we assumed a normal prior distribution cen- 
tered on zero for each of them. Practically, we added the follow- 
ing Bayesian penalty to our merit function: 



BPu 



- y (nXi\ 



(1) 



where TTVi is the TTV for the i' h C0R0T transit, and o-jt, is the 
error on its timing estimated by a preliminary individual analysis 
of this transit. 



5.2. Results 

Table 5 present the CoRoT-12 system parameters and l-cr error 
limits derived from our nominal MCMC run {MCMC\ ), and for 
the two other MCMC runs. 

Our MCMC analysis presents CoRoT-12b as an inflated 
Jupiter-mass planet (M p = 0.92 ± 0.07 M Jup , R p = 1.44 + 0.13 
Rj up ) transiting a solar analog star (M t - 1.08 + 0.08 M , 
R* - 1.1 ±0.1 R Q ). Using the stellar density deduced from our 
MCMC analysis (p, - 0.77 +JJ 2 " p ) and the effective temper- 
ature and metallicity obtained from spectroscopy (Table 3), a 
stellar evolution modeling based on the code CLES (Scuflaire et 
al. 2008) leaded to a stellar mass of 1.07 ± 0.10 M , in excellent 
agreement with our MCMC result, and to a poorly constrained 
age of 6.3 + 3.1 Gyr. It is also worth noticing that the two inde- 
pendent values obtained for the stellar surface gravity from our 
spectroscopic and global analysis are in good agreement (1.4 cr), 
indicating the good coherence of our final solution. 

Fig. 5 presents the period-folded C0R0T photometry binned 
per two minutes time intervals with the best fit transit model su- 
perimposed. The standard deviation of the residuals of this lat- 
ter LC is 592 ppm, demonstrating the excellent quality of the 
C0R0T photometry. 

Our results show that the limb-darkening coefficients u\ and 
u 2 are poorly constrained by the C0R0T transit photometry, de- 
spite its good precision. Indeed, the posterior distributions of u\ 
and u 2 are close to the prior distributions, indicating that the 
C0R0T data do not constrain these parameters much. 

Our final precisions on the stellar and planetary masses and 
radii are not excellent (about 7% on the masses and about 10% 



2 Jump parameters are the parameters that are randomly perturbed at 
each step of the MCMC. 



3 The stellar calibration law presented by Torres et al. is in fact func- 
tion of T e ff, [Fe/H] and logg. We modified it to use as input the stellar 
density instead of the stellar surface gravity (see Anderson et al. 2010b). 
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MCMC, 



Data 



Jump parameters 



MCMCi CoRoT transits 

HARPS (EGGS +H AM) 
HIRES (not transit) 



(R„/R s ) 2 , W, V, 
P, T ,K 2 ,Ci,C2, 

e cos a>, e sin a> 



MCMC 2 idem MCMCi 

+ CoRoT occultations 



idem MCMCi 
+ occultation depth dF 2 



idem MCMCi 



idem MCMCi 
+ 47 TTVs 



Normal prior distributions 



ill ~ N(0A7,0.03 2 ) 
u 2 ~ N(0.22, 0.02 2 ) 



m ~ N(0.47,0.03 2 ) 
u 2 ~ N(0.22, 0.02 2 ) 



m ~N(0.47,0.03 2 ) 

w 2 ~ N(0.22,0.02 2 ) 

TTV m - A1] ~N(0,aj, Tii ) 



Table 4. Specificities of the three MCMC runs performed during our global analysis. See text for details 



on the radii), and more observations are required to thoroughly 
characterize this system. In this context, improving significantly 
the precision on the stellar density (about 20%) is desirable. 
Such an improvement could be achieved mostly through a bet- 
ter characterization of the orbital parameters e cos a> and e sin to 
with more RV measurements (and possibly occultation photom- 
etry). Indeed, an new MCMC analysis assuming a perfectly cir- 
cular orbit leads > 2 times smaller error bars on the planet's and 
star's radii. The characterization of the system would also bene- 
fit from an improved determination of the transit parameters with 
more high-precision transit photometry, if possible acquired in a 
redder bandpass (less significant limb-darkening). 

The results of the run MCMC2 show that the occultation of 
the planet is not detected in the CoRoT data. We can only put an 
upper limit on its depth (3-cr upper limit = 680 ppm). 

As expected, the errors on To and P are significantly larger 
for the run MCMC3, but the posterior distributions obtained for 
the other parameters agree well with ones of the other MCMC 
runs. The resulting TTVs are shown in Fig. 6. No transit shows a 
significant timing variation. Still, the resulting TTV series seems 
to show a correlated structure. Fitting a sinusoidal function in 
this series leads to a best-fit period of about 24 epochs , i.e. of 
about 68 days. Nevertheless, the resulting false alarm probability 
is high, about 15%, indicating that this correlated structure is 
not very significant. Still, it is interesting to notice that, if we 
assume a rotational period of 68 days for the star and sin i - 
1, and using R, - 1.1 R e , we obtain a value of 1.2 kms^ 1 for 
v sin i, in excellent agreement with the value derived from our 
spectroscopic analysis (see Table 3). In this context, a possible 
interpretation of the low-amplitude structure visible in the TTV 
series is that it is caused by the rotation of the surface of the star 
and its influence on the transit barycenters. 

6. Discussion 

The position of CoRoT-12b in a planetary mass-radius diagram 
is shown in Fig. 7. While being denser than the extremely 
inflated planets WASP-17b (Anderson et al. 2010a), TrES-4b 
(Mandushev et al. 2007) and WASP- 12b (Hebb et al. 2009), 
CoRoT-12b appears to be a very low-density 'hot Jupiter'. Using 
the hypothesis that the planet is a core-less gazeous planet of so- 
lar composition, we used the planetary evolution code CEPAN 
(Guillot & Morel, 1995) to assess the ability of standard irra- 
diated planet models to explain the low-density of CoRoT-12b. 
Several models were used: a standard model with no extra heat 
source, a model for which the opacities were artificially mul- 
tiplied by 30, and three models with a constant energy deposit 
(10 26 , 10 27 and 10 28 erg.s -1 ) at the planet's center. Our results in 





dT [days] 

Fig. 5. Top: CoRoT transit photometry period-folded and binned 
per 2 minutes time intervals, with the best fit transit model su- 
perimposed. Bottom: residuals. Their standard deviation is 592 
ppm. 



terms of planetary size evolutions are shown in Fig. 8. For recall, 
we constrain the age of the system to 6.3 + 3.1 Gyr. Considering 
this age, the measured size of CoRoT-12b is in good agreement 
with all four evolution models. At most can we notice that an ex- 
tra heat source and/or of larger opacities are favored by the data, 
but a more precise radius measurement is needed to conclude. 

In this context, it is worth noticing that the precision on the 
planet's radius is mostly limited by the precision on the orbital 
eccentricity and argument of periastron (see Sec. 5.2). It is thus 
desirable to obtain more RV measurements of the system. Better 
constraining the planet's orbital eccentricity would also make 
possible the assessment of its past tidal evolution and its influ- 
ence on its energy budget (e.g., Ibgui et al. 2010). 
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Fig. 6. Top: median value and 1-cr limits of the TTV posterior 
distributions obtained in MCMC3. Bottom: same curve obtained 
after binning the TTVs per three (error of each bin = error on the 
mean). 
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planets (black circles, values from http://exoplanet.eu) in a mass- 
radius diagram. The error bars are shown only for CoRoT-12b 
(CI 2), WASP- 17b (W17), TrES-4b (T4), and WASP- 12b (W12) 
for the sake of clarity. 
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Fig. 8. Evolution of the size of CoRoT-12b (in Jupiter units, 
Rj U p) as a function of age (in billion years), compared to con- 
straints inferred from C0R0T photometry, spectroscopy, radial 
velocimetry and CESAM stellar evolution models. Red, blue 
and yellow-green dots correspond to the planetary radii and ages 
that result from stellar evolution models matching the inferred 
p*-T e s uncertainty ellipse within lcr, 2<x and 3<x, respectively. 
Planetary evolution models for a planet with a solar-composition 
envelope and no core are shown as plain lines and are labeled as 
follow : standard (black): irradiated planet with no extra heat 
source ; opacities x 30 (blue): opacities have been artificially 
multiplied by 30 compared to standard model ; dissipation: mod- 
els in which 10 26 (green), 10 27 (pink), and 10 28 (purple) ergs -1 
is deposited at the planet's center. These models assume a total 
mass of 0.92 Mj up and an equilibrium temperature of 1450 K 
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Parameter 



MCMCi 



MCMC 2 



MCMd 



Jump parameters 



Planet/star area ratio (R p /R s ) 2 

V = a cos i/R, [R t ] 

Transit width W [d] 

T - 2450000 [HJD] 

Orbital period P [d] 

RVA: 2 [ms-'d" J ] 

ecos w 

esinw 

Cl 
C2 

dF 2 



0.01744! 
0.609+1 



0.10726 
4398.62707 


+0.00089 
-0.00090 

± 0.00036 


2.828042 ± 


0.000013 


177!| 2 
-0 012+ 0024 

u,u -0.028 

053+°'° 73 

-0.064 

1.152 ±0.056 


0.028 ± 


0.052 



n nnW'-" 0044 

u.ui;j7_ 000041 

614+ 0060 

"• ui -0.056 

1071+ 00013 

u.iu; i_ ooon 

4398 62704+°°°° 38 
9 sosnzr?+oooooi3 

/ -° /8U4J -0.000014 

176 ±11 
000+° ° 2 ° 

"•"""-0.040 

069+° ° 69 

u,uu -0.082 
i i co+0.054 
1,1JJ -0.059 

031+ 0052 

"■"■"-0.051 

00009+° °° 022 

U.UUW7J nnnno 



U.U1 '-->- , _o.00043 

n SQ9+ 0040 
"• J5 ^-0.046 

1071 +00 °" 
4398 6266+ 00013 

9 STSflft r+0.000052 

177 ± 10 
-0 017+ 0024 

"•"^ -0.026 

043+ 0072 

1 146 +0.058 

027+ - 048 

"•"^ -0.049 



Deduced stellar parameters 



U\ 

Density p, [p Q ] 

Surface gravity log g, [cgs] 

Mass M, [M e ] 

Radius R, [R Q ] 



0.466 ± 0.027 

0.219 ±0.021 

77+ " 20 

"•' '-0.15 
a 175+0.065 
J -0.062 

1 ' U/8 -0.072 

1 116+ 00% 

1,llu -0.092 



0.468!°™ 

0.217 ± 0.020 

75+ a2 ° 

"■' J -0.15 

4.366+°°g 
1.083+°;™ 

1 17Q+0.097 
1 - lz - 7 _0.092 



464+ 0028 

u.tut_ 0025 

n 9 i n+0.020 
"• z - 1? -0.019 

0.81!°- 

4-388!°;°S 

1 076+ 0077 

1# " '"-0.071 

1 098+ 0072 

1.VJ70 0Q76 



Deduced planet parameters 



RV «[ms-'] 

U 'transit l/^*J 
"occuhation L***J 

T occultatio „ - 2450000 [HJD] 

Orbital semi-major axis a [AU] 

Orbital inclination i [deg] 

Orbital eccentricity e 

Argument of periastron uj [deg] 

Equilibrium temperature T eq [K]" 

Density p p [p ]up \ 

Surface gravity log g p [cgs] 

Mass M p [M Jup ] 

Radius R p [R Jup ] 
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573+ 0027 
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0.04016" 
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1 m i +0.083 
-'• v " i -0.077 
Q 916 +0.068 
"• 7iu -0.064 

1 45+ 013 



3 060+ 0065 

f) Q 1^+0.068 
"■ 71J -0.064 

1 41+0.10 
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Table 5. Median and 1 -cr limits of the posterior distributions obtained for the CoRoT- 1 2 system derived from our three MCMC runs 
(see Table 4). MCMC\ is our nominal analysis (see text for details). "Assuming A=0 and F—\. 



